High-field quantum calculation reveals time-dependent negative Kerr contribution 
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The exact quantum time-dependent optical response of hydrogen under strong field near infrared 
excitation is investigated and compared to the perturbative model widely used for describing the 
effective atomic polarization induced by intense laser fields. By solving the full 3D time-dependent 
Schrodinger equation, we exhibit a supplementary, quasi- instantaneous defocusing contribution miss- 
ing in the weak-field model of polarization. We show that this effect is far from being negligible in 
particular when closures of ionization channels occur and stems from the interaction of electrons 
with their parent ions, ft provides an interpretation to higher-order Kerr effect recently observed in 
various gases. 

PACS numbers: 34.80.Dp, 42.65.An 



Atoms and molecules exposed to a nonresonant intense 
(10-100 TW/cm^) laser field exhibit highly nonlinear dy- 
namics that has motivated a wealth of experimental and 
theoretical studies and led to the observation of phenom- 
ena such as above-threshold ionization (ATI) [Jj, high- 
harmonics generation (HHG) attosecond pulses gen- 
eration 0, or filamentation [JQ- Understanding the 
first three processes has required to describe the atomic 
and molecular dynamics under strong-field excitation. In 
this regime, the electric field is as intense as the intra- 
atomic field, so that electronic transitions can be en- 
hanced by dynamic resonances leaving the atom in a co- 
herent superposition of bound and continuum states. In 
the same time, ionization channel closures occur, i.e. the 
minimum number of photons required to ionize the atom 
increases by several units due to the Stark shift of the ion- 
ization potential Because of this complex dynamics, 
the atomic optical response can no longer be described 
as a perturbative series of the field. Instead, the atom 
dynamics and its associated polarization must be evalu- 
ated at each time during the interaction by solving the 
time-dependent Schrodinger equation (TDSE) describing 
the interacting atom. As far as the filamentation process 
is concerned, it is generally interpreted through a pertur- 
bative approach as resulting from a dynamic balance be- 
tween Kerr self-focusing and defocusing by the free elec- 
trons originating from ionization during the interaction 
(Drude model). Retrospectively, it is surprising that this 
process is still described in the perturbative framework, 
contrary to HHG or ATI. Indeed, since the typical in- 
tensity in filaments (50 TW/cm^ 0, ^) is of the same 
order as that used in HHG or ATI experiments, it seems 
natural to wonder whether filamentation can still be ac- 
curately described within the lowest order perturbation 
theory, i.e. using the commonly used Kerr and Drude 
(KD) model. In particular, the saturation and even the 
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inversion of the Kerr effect (i.e. the sign inversion of 
the nonlinear refractive index), empirically described as 
negative higher-order Kerr effect (HOKE) , were reported 
between 19 and 33 TW/cm^ in major air components 
[Til , [l^ . A laser- induced transient grating [l^ was also 
proposed as an alternative interpretation of the measure- 
ments, casting doubts on the relevance of this model. 
Nevertheless, saturation and inversion of the Kerr index, 
predominantly attributed to continuum-continuum tran- 
sitions, were predicted in hydrogen [343 j in atomic 
silver [l3, in Hj [l3|, as well as with TDSE calculations 
in a one-dimensional Dirac potential [l9| , or by consider- 
ing Kramers-Kronig relations and multiphoton ionization 
rates 

In this Letter, we perform 3D non-perturbative ab- 
initio calculations of the interaction between a near in- 
frared strong laser field and an hydrogen atom. The 
choice of hydrogen is motivated by the possibility of per- 
forming rigorous exact calculations without any assump- 
tion, in particular about the exact form of the potential. 
We show that the usual perturbative model of polariz- 
ability, describing a purely third-order Kerr medium to- 
gether with a negative plasma contribution following the 
Drude model, misses a substantial additional negative 
contribution associated with ionization channel closures. 
The latter contribution is due to the over-acceleration 
of the electrons interacting with their parent ions dur- 
ing the pulse and can be identified as the empirically- 
introduced higher order Kerr effect (ill . [T2| . This effect 
may affect any strong-field experiment in long gaseous 
media. The Kerr effect is further influenced by field- 
induced resonances that eventually leave several excited 
states significantly populated at the end of the pulse. 
Moreover, we show that the Kerr inversion mechanism 
cannot be specifically attributed to neither bound-bound 
nor continuum-continuum transitions that are gauge de- 
pendent. This stems from the fact that the atomic states 
are not eigenstates of the dressed atomic system, and 
consequently are not quantum observables of the system 
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during the interaction. This result sheds a new hght on 
the higher order Kerr effect controversy which liad crys- 
tahzed the debate around this question (f^ - [l8| . 

Within the dipole approximation, the TDSE describing 
the electron wavefunction evolution in the presence 
of an electric field E(i) linearly polarized along the axis 
z reads: 



\g=23 TW/cn? 



Ijj=32.5 TW/cm^ 



(1) 



where Hq = — f/r is the hydrogen atom Hamil- 

tonian, iJint ~ ^E(t) • r (resp. Hint = -A-(t) • tt, where 
A{t) is the vector potential such that E{t) = —dA/dt 
and TT = — iV) is the interaction term expressed in the 
length gauge (resp. velocity gauge). The time-dependent 
wavefunction \ip) is expanded on a finite basis of B-splines 
PH allowing memory efficient fast numerical calculations 
with a very large basis set: 
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Figure 1. (Color online) Nonlinear polarization envelope as 
a function of time for a 30 (a), 32.5 (b), 38 (c), and 42 (d) 
TW/cm^ 48 cycles pulse. 
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where and Y"'' are B-spline functions and spheri- 
cal harmonics, respectively. The basis parameters (Imax, 
'T-max, k and the spatial box size) and the pro pag ation pa- 
rameters are chosen to ensure convergence [23 ■ Unless 
otherwise specified, calculations are performed in the ve- 
locity gauge where computation is eased The atom 
is initially in the ground state (Is) and the electric field E 
is expressed as E{t) = cos(t/(Tt)^ cos(cjoi), where ojq is 
the central pulsation of the laser, at = 2N/ujo, and N the 
total number of optical cycles in the temporal window. 
The simulations are performed for a laser wavelength of 
800 nm and pulse durations (FWHM of intensity) of 23, 
29, 47, 58, and 93 fs (corresponding to iV=24, 30, 48, 60, 
and 96 cycles respectively). 

The microscopic polarization p(t) and the medium po- 
larizability a at the driving frequency luq are defined in 
atomic units as: 



Pit) 
a{uJo) 



{if; {t)\r cos 9\il;{t)) 



(3) 
(4) 



where p(w) (resp. E{lu)) is the Fourier transform of p(t) 
(resp. E{t)) and 9 is the angle between the dipole mo- 
ment and its projection along the electric field. No collec- 
tive effect is expected at the considered time scale (23| . 
even if a comprehensive treatment would need to con- 
sider the interaction with the potential of the neighbor- 
ing atoms, which is beyond the scope of the present work. 
Under the isolated atom assumption, the macroscopic po- 
larization P is then calculated as P{t) = Afp{t), with TV 
the number density. 



Its nonlinear part Pnl and the nonlinear polarizability 
Aa are calculated as: 



PNL(t) = P{t)-Af lim a{Io)E{t), (5) 
Aa{Io) = a{Io) - lim a(/o), (6) 

In M-O 



where Iq is the peak intensity of the pulse. As shown in 
Fig. [l] the nonlinear polarization envelope (the carrier 
and the harmonics have been filtered out) follows the in- 
tensity profile in the low-field regime (/o<30 TW/cm^), 
as expected for a purely cubic Kerr effect. As the inten- 
sity increases, a transient, time-retarded negative contri- 
bution develops on the second half of the pulse and leads 
to a sign inversion of the nonlinear polarization in the 
falling edge of the pulse. In the frequency domain, the 
nonlinear polarizatiblity of the atom at the fundamental 
frequency luq increases linearly with the pulse peak inten- 
sity, saturates, and eventually becomes negative (Fig. (2). 
The threshold intensity for inversion is in line with previ- 
ous works [Hill , [l^ and similar for all investigated 
pulse durations. 

From this point, one can wonder if the usual way to 
describe the atomic polarizability at the fundamental fre- 
quency through a perturbative approach still remains 
valid at high intensity or if supplementary effects must 
be considered. To answer we shall first recall that the 
KD model of polarization is based on the following as- 
sumptions: i) atoms are responsible for a pure cubic Kerr 
effect, ii) free electrons are produced through ionization 
during the interaction with the field while depletion of the 
ground state population remains negligible, iii) the ion- 
ized electrons (commonly called plasma), insensitive to 
the atomic potential, accelerate in the field as free par- 
ticles and are responsible for balancing the Kerr effect. 
Following the above assumptions, the effective nonlinear 
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Figure 2. (Color online) Nonlinear polarizability as a function 
of peak intensity for 24 (a), 48 (b), 60 (c) and 96 (d) cycles 
pulse calculated with TDSE (red solid line) and KD model 
(blue dashed line), respectively. 



macroscopic polarization of the system reads: 



oeff 
NL,KD 



Pc 



E{t), 



(7) 



where p is the free electrons density, pc = nicUJQeo/q^ 
is the critical plasma density, n2 is the intensity inde- 
pendent nonlinear index of the medium (which can be 
extracted from TDSE calculations at weak intensity), eo 
is the vacuum permittivity, and mo and q are the elec- 
tron mass and charge, respectively. Since Eq. [7] depends 
on the free electrons density, the perturbative scenario 
de facto assumes that the latter is a quantum observable 
during the interaction, and, as such, can be measured 
and evaluated. In that framework, models like multipho- 
ton ionization, PPT, or ADK can be used to evaluate the 
ionization probability during the interaction (23 - l26| . It 
should be emphasized that this approach is inappropri- 
ate when dealing with quantum mechanics because the 
wavefunction describing the system depends on the rep- 
resentation of the interaction (i.e. the gauge), whereas 
the physical predictions (for instance, the electron po- 
sition and mechanical momentum) do not. As shown 
in Fig. [3Ja), the probability for the electron to be in 
the continuum of the unperturbed atom (correspond- 
ing to the ionized electrons density) varies by several 
decades during the interaction depending on the consid- 
ered gauge and therefore is not a quantum observable. 
The ionized electron density cannot therefore be mea- 
sured during the interaction, making the time-dependent 
KD scenario questionnable from a quantum mechani- 
cal point of view. To show the divergence between the 
full TDSE calculations and the KD model, we used the 
Perelomov- Popov- Terent'ev (PPT) model [23 for evalu- 
ating p. Since PPT overestimates the postpulse electron 
densities calculated with TDSE, which gives exact results 
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Figure 3. (Color online) (a) Probability of the electron to be in 
the continuum during a 48 cycles 10 TW/cm^ pulse calculated 
in the length and velocity gauges, respectively, compared with 
the output of the PPT formula, (b) Correction factor K of 
the PPT ionization probability. Nonlinear polarizability as a 
function of intensity calculated in (c) the length and (d) ve- 
locity gauge, respectively. Aobb, Aqbc, and Aacc refer to 
the partial nonlinear polarizabilities related to bound-bound, 
bound-continuum, and continuum-continuum transitions, re- 
spectively. 



in the case of atomic hydrogen, we corrected it so that 

p{t = +Oo) = PTDSE(t = +Oo) = Kp-pvTit = +Oo), 

where K is an intensity dependent correction factor de- 
picted in Fig. El^b). Note that, the amount of ionized 
electrons becomes an appropriate quantum observable af- 
ter extinction of the external electric field [see Fig. ISJa)]. 

The KD model qualitatively reproduces the nonlinear 
polarization dynamics calculated with TDSE (Figs. [T] 
and [2]). However, it significantly under-estimates the 
negative contribution leading to the sign inversion of the 
nonlinear polarizability. As shown in Fig. [21 the discrep- 
ancy between the two models develops as soon as the 
intensity exceeds the 10-photons ionization channel clo- 
sure (/q ~ 32 TW/cm^. see Fig. [2]) where 11 photons are 
necessary to ionize the atom as the ionization threshold 
shifts upward due to the AC stark effect, and further in- 
creases above the 11-photons ionization channel closure 
at Iq =57 TW/cm^. Ionization channel closure turns out 
to be a key ingredient of the negative nonlinear transient 
contribution to the polarization. Its impact appears sim- 
ilar to that induced by the ionization suppression mech- 
anism (27| but takes place at lower intensities. 

During the interaction, the KD model of polarization 
deviates from the quantum calculations with the fourth 
power of the incident peak intensity, consistent with its 
empirical identification with a negative ng/g contribu- 
tion to the refractive index [Tl] | at large intensity. The 
negative effect induced by ionization channel closures is 
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Figure 4. (Color online) Spectral dependence of the nonlinear 
polarizability for (a) 24 and (b) 48 cycles. 



therefore similar (except for the delay of few femtosec- 
onds, compatible with the experimental temporal resolu- 
tion), to the higher order Kerr effect introduced empiri- 
cally to describe the experimental observations of intra- 
pulse Kerr saturation and inversion [Til. IT^. The major 
air constituents, as well as most usual gases, exhibit ion- 
ization channel closures in the same intensity range as hy- 
drogen (21.5 and 48 TW/cm^ in Ar, 32 and 57 TW/cm^ 
in O2, and 25 and 51 TW/cm^ in N2). As a consequence, 
in all of these gases, the effect of ionization channel clo- 
sures has to be taken into account in the sign reversal of 
the nonlinear refractive index. It should be noted that 
the nonlinear polarizability saturates and reverses its sign 
only at frequencies within the bandwidth of the incident 
pulse spectrum (see Fig. S]) as recently observed in (28j . 
This unexpected finding is consistent with the apparent 
contradiction between studies performed at a single wave- 
length, which exhibit refractive index saturation and in- 
version fill . [l^. 1 2 81] ■ and two-color investigations, that do 
not [H-lsoi. 

Beyond the ionization channel closure, excited states 
dynamically shift with intensity. The resulting dynami- 
cal resonances lead to a net population transfer to some 
excited states after the pulse turn-off. For instance, 2d, 
3d, and l-4g states are signifiquantly populated after the 
interaction with a pulse of ~40 TW/cm^ peak intensity, 
impacting the atomic polarization [Fig. [Ud)] and the as- 
sociated nonlinear refractive index. This time-dependent 
effect becomes even more pronounced for longer pulses 
(>50 fs) where the system remains resonant over long 
time. It may therefore seem natural to investigate 
which transitions (e.g., bound-bound, bound-continuum, 
and continuum-continuum transitions as defined in (l^ ) 
mostly contribute to the Kerr saturation and inversion 
in order to exhibit a simple physical interpretation about 
the sign inversion of the Kerr effect at high intensity. 
However, like the amount of ionized electrons in the 
presence of the external field, the bound-bound, bound- 
continuum, and continuum-continuum contributions to 
the total polarizability (i.e. the partial polarizabilities) 
are also gauge-dependent, as illustrated in Figs. [3jc,d). 
This prevents any rigorous identification of specific tran- 
sitions responsible for the Kerr saturation and inver- 
sion |l4l4l7l |. More particularly, while only continuum- 



continuum transitions (i.e. the acceleration of electron in 
the continuum) seem to be responsible for the Kerr sign 
inversion in the length gauge, both continuum-continuum 
and bound-bound transitions actively participate in the 
velocity gauge. The fact that the partial nonlinear polar- 
izability Aacc induced by continuum-continuum transi- 
tions is a major negative contribution to the total non- 
linear polarizability in both velocity and length gauges 
does not imply any similar behavior in other gauges. As 
a consequence, while the higher order Kerr effect contro- 
versy has crystalized the debate on whether or not ion- 
ized electrons (or equivalently the plasma) are the main 
source of the Kerr saturation and inversion mechanism, 
it appears that this question is in fact irrelevant due to 
the gauge-dependency of the answer. 

As a conclusion, time dependent non perturbative cal- 
culations show that the nonlinear polarizability of hydro- 
gen saturates and reverses its sign around 30 TW/cm^, 
an intensity at which refractive index inversion is ob- 
served 111, l2|. In this strong-field regime, characterized 



by ionization channel closure, the electrons promoted 
into the continuum significantly interact with their par- 
ent ions and cannot be treated as free charged particles 
during the pulse. As a consequence, the well admitted 
scenario of Kerr effect saturation and inversion misses a 
quasi-instantaneous supplementary negative term, which 
can be identified as higher order Kerr effect or alterna- 
tively as a deviation from the Drude model. This result is 
compatible with the original observations of Kerr satura- 
tion and inversion (ill . [l3 | , even if plasma grating effects 
could have a partial contribution in the experimental ob- 
servations. It is worth mentioning that a direct compar- 
ison with the birefringence experiments (ill , [l^ would 
call for an evaluation of the refractive index along two 
orthogonal directions. Such a calculation would need to 
include an electric field polarized along these two direc- 
tions as well as magnetic sub levels transitions, which is 
beyond the scope of this work. 

Moreover, seeking to pinpoint specific transitions be- 
tween atomic states (either bound or in the continuum) 
as responsible for the nonlinear refractive index inversion 
is irrelevant, since the intra-pulse populations (in partic- 
ular, the free electron density or plasma density) are not 
physical observables, but rather gauge-dependent quan- 
tities. Further studies will be needed to exhibit a simple 
and gauge independent physical interpretation of the sign 
inversion of the Kerr effect. These results therefore draw 
a radically new perpective on the current debate about 
higher order Kerr effect. 

Although derived in atomic hydrogen, these conclu- 
sions should apply to all major air constituents since 
ionization channel closures occur in the same intensity 
regime in all of these gases. Further works will be needed 
to describe the Kerr effect with a ready-to-use, computa- 
tionnally efficient parametrization. In particular, prelim- 
inary studies suggest that the full quantum trajectory 
of the electron could be described as those of a parti- 
cle with a pulse intensity- and shape-dependent effective 
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mass. This parametrization will be useful for practical 
purposes requiring repeated computation of the nonlin- 
ear refractive index, like the simulation of laser filamen- 
tation. 

Acknowledgments. This work was supported by the 



[1] J.H. Eberly, J. Javanainen, and K. Rzazewski, Phys. 

Rep. 204 (5), 331-383 (1991) 
[2] M. Ferray et al, J. Phys. B 21, L31 (1988) 
[3] P. Antoine, A. L'Huillier, and M. Lewenstein, Phys. Rev. 

Lett. 77 (7), 1234-1237 (1996) 
[4] S. L. Chin et al, Can. J. Phys. 83, 863 (2005). 
[5] L. Berge et al., Rep. Prog. Phys. 70, 1633 (2007). 
[6] A. Couairon and A. Mysyrowicz, Phys. Rep. 441 47 

(2007) . 

[7] J. Kasparian and J. -P. Wolf, Opt. Express 16, 466 

(2008) . 

[8] R. Wiehle, B. Witzel, H. Helm, and E. Cormier, Phys. 

Rev. A 67, 063405 (2003) 
[9] J. Kasparian, R. Sauerbrey, S. L. Chin, Applied Physics 

B 71, 877 (2000) 
[10] A. Becker et a/.. Applied Physics B 73, 287 (2001) 
[11] V. Loriot et al., Opt. Express 16, 13429 (2009); Erratum 

in Opt. Express 18 3011 (2010) 
[12] V. Loriot et al., Laser Physics 21 (7), 1319-1328 (2011) 
113] J.K. Wahlstrand and H.M. Milchberg, opt. Lett. 36(19), 

3822-3824 (2011) 
[14] M. Nurhuda, A. Suda, and K. Midorikawa, Phys. Rev. 

A 66, 041802 (2002) 
]15] M. Nurhuda, A. Suda, and K. Midorikawa, New J. Phys. 



Conseil Regional de Bourgogne (FABER program), the 
FASTQUAST ITN Program of the 7th FP. JPW ac- 
knowledges financial support from the ERC advanced 
grant "FilAtmo". EC and PB thank O. Peyrusse and 
the CRI-CCUB for CPU loan on their respective multi- 
processor servers. 



10, 053006 (2008) 
]16] P. Kano, M. Brio, and J.V. Moloney, Comm. Math. Sci. 

4 (1), 53-80 (2006) 
]17] E.A. Volkova, A.M. Popov, and O.V. Tikhonova, JETP 

Lett. 94 (7), 519-524 (2011) 
]18] E. Lorin et al, Physica D 241(12), 1059-1071 (2012) 
[19] A. Teleki, E.M. Wright, and M. Kolesik, Phys. Rev. A 

82, 065801 (2010) 
[20] C. Bree, A. Demircan, and G. Steinmeyer, Phys. Rev. 

Lett. 106, 183902 (2011) 
[21] H. Bachau et al. Rep. Prog. Phys. 64, 1815-1942 (2001) 
[22] E. Cormier and P. Lambropoulos, J. Phys. B 29(9), 1667- 

1680 (1996) 

[23] B. Pasenow et al., Opt. Express, 20, 2310-2318 (2012) 
[24] L. V. Keldysh, Sov. Phys. JETP 20(5), 1307-1314 (1965) 
125] A. M. Perelomov, V. S. Popov, and M. V. Terentev, Sov. 

Phys. JETP 23(5), 924-934 (1966) 
126] M. V. Ammosov, N.B. Delone, and V. P. Krainov, Sov. 

Phys. JETP 64(6), 1191-1194 (1986) 
127] F. Morales et al, PNAS 108, 16906 (2011) 
128] J. H. Odhner et al., Phys. Rev. Lett., 109, 065003 (2012) 
129] J. K. Wahlstrand et al., Phys. Rev. Lett. 107, 103901 

(2011) 

130] J. M. Brown et a/., Opt. Lett. 37, 1064 (2012) 



